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Abstract 

We consider the mixed problem for L the Lame system of elasticity in a 
bounded Lipschitz domain f2 c R 2 . We suppose that the boundary is written 
as the union of two disjoint sets, d£l = D U N. We take traction data from 
the space L P (N) and Dirichlet data from a Sobolev space W 1,P (D) and look 
for a solution u of Lu = with the given boundary conditions. We give a scale 
invariant condition on D and find an exponent po > 1 so that for 1 < p < po, we 
have a unique solution of this boundary value problem with the non-tangential 
maximal function of the gradient of the solution in L p (dQ). We also establish 
the existence of a unique solution when the data is taken from Hardy spaces 
and Hardy-Sobolev spaces with p in (pi, 1] for some p\ < 1. 



1 Introduction 

We consider the L p -mixed problem for the Lame system of elasticity in a Lip- 
schitz domain Q, C R 2 . Thus we consider the operator L = + (A + /i)V div 
acting on vector- valued functions. We assume that the Lame parameters sat- 
isfy // > and A > —fi so that the operator L is elliptic. (See the discussion 
after ()2.5[) for more details.) We assume that we have a decomposition of the 
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boundary into two sets d£l = D U N where D C dtt is a non-empty, open, 
proper subset of d£l and ./V = dfl\D. We let du/dp denote a traction operator 
at the boundary. Given data /jv on iV and fjj on D, we look for a solution u 
of the boundary value problem 



Here, (Vn)* denotes the non-tangential maximal function of Vtt. The goal of 
this paper is to give conditions on D and the data fu and that will allow us 
to establish the existence and uniqueness of the solution u and show that the 
solution u depends continuously on the data. 

This investigation is a continuation of work that dates back at least to 
Dahlberg [8] who studied the Dirichlet problem for the Laplacian in a Lipschitz 
domain by a careful investigation of harmonic measure. Further developments 
for the Laplacian include Jerison and Kenig [19] who treated the Neumann 
problem and the regularity problem in L 2 (dft) and Dahlberg and Kenig [9] 
who studied the Neumann and regularity problems in L p {d£l) for p between 
1 and 2. We point out that it has become standard to refer to the Dirichlet 
problem with data that has one derivative in some L p (dQ) space as the regu- 
larity problem. In both the regularity and the Neumann problem, the goal is 
to establish estimates on the non-tangential maximal function of the gradient 
of the solution. Other relevant developments include the work of Dahlberg, 
Kenig, and Verchota [11] who studied the Dirichlet, Neumann, and regularity 
problems for the Lame system in L 2 (dQ) and work of Dahlberg and Kenig 
|10j who studied the regularity and Neumann problems for the Lame system in 
L v (dQ) in dimension 3. Additional results were obtained for the Lame system 
by Mendez and M. Mitrea [22] and Mayboroda and M. Mitrea [21] and again 
their work is limited to low dimensions. 

Another strand of this story is the study of the mixed problem in Lipschitz 
domains. The author Brown and collaborators including Capogna, Lanzani, 
and Sykes [2[ [20], [29] [30] have established well-posedness of the mixed prob- 
lem for the Laplacian in restricted classes of Lipschitz domains. I. Mitrea 
and M. Mitrea [23] have extended these results to cover the Poisson prob- 
lem for the Laplacian in a wide range of function spaces. These methods 
have been extended to the Lame system with I. Mitrea [4], the Stokes sys- 
tem with I. Mitrea, M. Mitrea, and Wright [I], and to the Hodge Laplacian by 
Gol'dshtein, I. Mitrea, and M. Mitrea [18]. The aforementioned works rely on 
the use of the Rellich identity with a vector field a chose that a ■ v changes 
signs as we pass from D to N . Here v is the outer unit normal to d£l. It 
seems unlikely that this technique can be applied in a general Lipschitz do- 
main. In two recent papers, the authors and Taylor [26] [33] have developed 
techniques to investigate the mixed problem for the Laplacian in a general Lip- 
schitz domain and for quite general decompositions of the boundary. Another 
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interesting approach was taken by Venouziou and Verchota [35] who study the 
mixed problem for the Laplacian in a class of polyhedral domains. The Ph.D. 
dissertation of Venouziou [M] includes a discussion of mixed problems for the 
biharmonic equation in a class of polyhedral domains and some examples where 
these mixed problems are not solvable. Venouziou also gives a nice summary 
of early work on the mixed problem. 

The work reported here extends the methods of Taylor, Ott, and Brown to 
the Lame system in two dimensions. A key ingredient of our approach to the 
mixed problem (11.11) is estimates for the Green function for the mixed problem. 
Taylor, Ott, and Brown [314 [33] obtain the needed estimates for the Green 
function for the Laplacian in all dimensions by the method of de Giorgi [12J. (As 
we are considering constant coefficient operators, solutions are well-behaved in 
the interior of the domain. The interesting issue is regularity at the boundary.) 
However, this approach is not available for elliptic systems. The Green function 
for the mixed problem in two dimensions was studied recently by Taylor, Kim, 
and Brown |32j . We note that the work of Dahlberg and Kenig jlO] provides 
estimates for the Green function for the Dirichlet problem and the Neumann 
problem in three dimensions. However, this approach depends on estimates 
for the L^-Dirichlet problem for p near 2. It is not clear how to extend this 
approach to mixed boundary value problems. The investigation of the mixed 
problem for elliptic systems in higher dimensions remains an interesting open 
question. 

With the estimates for the Green function in hand, we argue as in the 
work of Dahlberg and Kenig [9] and use the asymptotic behavior of the Green 
function to study our boundary value problem for p < 1. This gives existence 
of solutions when the data comes from atomic Hardy spaces. Finally, we adapt 
an argument of Shen [28] to establish results in LP for p > 1 and obtain the 
following theorem which is the main result of this paper. 

Theorem 1.2 Let $7 C R 2 be a Lipschitz domain and let D C dfl be a non- 
empty open proper subset satisfying the corkscrew condition \2. 1\) and the con- 
ditions \2. 2\) and \2. 3\) . Let L be the Lame operator with coefficients satisfying 

(KM). 

1) There exists po > 1 so that for 1 < p < po, the LP -mixed problem is 
well-posed in the sense that if fN £ L P (N), fp € W 1,p (dQ,), the problem (COP 
has a unique solution which satisfies the estimate 

\\(y u )*\\Lp(dn) < C(II/aHIlp(ao + WfoWw^pidn))- 

The boundary values of u and Vu exist as non-tangential limits. 

2) There exists p\ < 1 so that if pi < p < 1, the LP -mixed problem is well- 
posed. Thus if fN £ LL P (N) and fu lies in the Hardy-Sobolev space H 1,p (d£l), 
then there exists a unique solution of the LP -mixed problem which satisfies 

du 

\\u\\H^p(dn) + ||-o-||ifi>(an) + ||(Vn)*||xj,(an) < C(\\f N \\ HP ( N ) + ||/D||ifi,P(an))- 
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We point out that we will require that our Dirichlet data have an extension 
from D to 30. The sets D that we consider may not be extension domains for 
Sobolev spaces. As our solutions (at least for p > 1) will have boundary values 
in W 1,p (dQ), it is clear that we need to assume that fo is the restriction to D 
of a function which lies in W l,p (d£l). 

Section [2] gives the definitions and assumptions needed in the rest of the 
paper. Section [3] gives a reverse Holder estimate for the gradient of a weak 
solution of the mixed problem. In section U we give the fundamental estimate 
for solutions with atomic data. Uniqueness is treated in section [5] and the 
proof of the uniqueness assertion of Theorem 11.21 is given in Theorem 15.81 The 
existence of solutions for the Hardy space problem with p < 1 is given in 
Theorem 16.11 of section [6] and the existence for p > 1 is indicated in Theorem 
16.31 which completes the proof of Theorem 1.2. Two appendices summarize 
information that is probably known, but not available in a convenient format. 
Appendix lAl gives several versions of the Sobolev inequalities that are essential 
to this work and Appendix |B] provides a treatment of the regularity problem 
with data from the Hardy-Sobolev space H ,p {dQ). 

2 Definitions 
2.1 Domains 

We will assume that is a bounded Lipschitz domain. While our final results 
hold only in two dimensions, many steps of the argument can be carried out 
in any dimension. When possible we will give our arguments in R n , n > 2, in 
order to highlight that part of the argument that is truly two-dimensional and 
to lay out the issues that restrict our argument to two dimensions. We assume 
O is a bounded open set in R n and for M > 0, r > 0, and x € 30, we define 
Z r {x) by Z r (x) = {y : \x' — y'\ < r, \x n — y ri \ < (4M + 2)r}. We say that Z r (x) 
is a coordinate cylinder for 30 if there is a Lipschitz function (j) : R n_1 — > R 
with Lipschitz constant M so that 



The coordinate system (y', y n ) £ R n_1 x R is assumed to be a rotation and 
translation of the standard coordinate system on R n . 

We say that is a Lipschitz domain with constant M if for each x € 00, we 
may find a coordinate cylinder for dVt centered at x, Z r (x). Since 30 is compact, 
we may find ro and a finite collection of coordinate cylinders {Z ro (xi)}^L 1 which 
cover 30 and so that for each i, Zioi ro (xi) is also a coordinate cylinder. 

Our estimates near the boundary will be defined using the non-tangential 
maximal function. To define this object, we fix a > and for x G 30, we 
define the non-tangential approach region with vertex at x by T(x) = {y G O : 
\x — y\ < (1 + a)d(y)} where d(y) = dist(y,<90) denotes the distance from y 



30 n z r (x) 
on z r (x) 



{y-y n = Hv')} n z r{x) 

{y:y n >4>{y')}C\Z r {x). 
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to the boundary. For a function v defined in fi, we define the non-tangential 
maximal function of v by 

v*(x) = sup \v(y)\, x G dd. 
y eT(x) 

It is well-known that the L p norm with respect to surface measure of the non- 
tangential maximal functions defined using different choices of a are compara- 
ble. Thus, we suppress the dependence of T and v* on the parameter a. 

We recall a useful tool from the Ph.D. dissertation of G. Verchota [36, 
Theorem A.l]. Given a Lipschitz domain C R n , we may find a family of 
domains {&k}kLi with Hj, C with the properties that a) each J)*, is a C°°- 
domain, b) there are bi-Lipschitz transformations : d£l — > d£lk such that 
Afc(x) converges to x. Using Verchota's construction we can see that the A^(x) 
converges to x non-tangentially as k — > oo, i.e. if the constant a used to define 
the approach regions T(x) is sufficiently large, then Aj~(x) G T(x) for all k. 

We define a star-shaped Lipschitz domain. These domains will be useful 
because they have sufficiently regular geometry that we can establish estimates 
that are uniform over the class of domains. Let x G R n , r > 0, M > 0, and 
4> : S n_1 — > [1, 1 + M]. We say that Q, is a star-shaped Lipschitz domain with 
center x, constant M, and scale r if (j> is Lipschitz with constant M and 

U = x + {y : \y\ < r<p(y/\y\)}. 

We introduce boundary balls (or intervals if n = 2). For x G d£l and 
r G (0, 100r ), we let A r (x) = Z r (x) flffi. If x G tt and r G (0, 100r ), we 
define local domains as follows. If d(x) = dist(x,9Q) > r, then Q r (x) = B r (x), 
the ball centered at x with radius r. If < r, then x lies in one of the 
coordinate cylinders Z ro (xi) = Z\. If x = (x',x n ) in the coordinate system for 
the cylinder Zj, then we put x = (x',4>(x')) and define Q r (x) = Z r (x). We let 
x* = (x',(j)(x') + (4M + l)re n ) and observe that the domains Q r (x) are star- 
shaped with respect to each point in B r / 2 {x*), thanks to our choice of M. This 
allows us to prove scale-invariant Sobolev inequalities and Korn inequalities in 
these domains. If x is on or near the boundary, then x may lie in several of 
the coordinate cylinders which cover the boundary and thus we have several 
choices for A r (x) and Q r (x). Our estimates will hold for any such choice with 
the condition that if several of these objects appear in an estimate we make 
a consistent choice of coordinate cylinder to define them. Note that in earlier 
works, we have sometimes used S r (x) = d£l n B r {x) and S r {x) = fl n B r (x), 
x G 0, in place of A r (x) and Cl r (x). Each set of definitions has its advantages 
and at one point in the argument below we will find it convenient to use S r (x) 
in place of A r (x). 

We specify our decomposition of the boundary by giving our assumptions 
on D. We assume D C d£l is a non-empty proper open subset of d£l. We will 
assume that D satisfies a corkscrew condition. To describe this condition, let 
A C dVt denote the boundary of D in dVt and 5{x) = dist(x, A). We say that D 
satisfies the corkscrew condition if for all x G A and r in (0, 100ro), there exists 
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x r £ D so that \x r — x\ < r and 5(x r ) > M _1 r. In Taylor, Ott, and Brown 
[331 Section 2], we show that if D satisfies the corkscrew condition, we have the 
following consequence: 

If x G D and r € (0, lOOro), then there exists x r € D so that , , 

\x — x r \ < r and A M -i r (x r ) C D . 

We note that the extreme cases D = or (90 of the mixed problem (II. ip 
correspond to the traction and regularity problems, respectively. As these have 
been studied elsewhere, we exclude these cases from our discussion in this paper. 

We also require that D satisfy the following conditions that involve the 
integrability of 5. To state these conditions, for r € (0, lOOro), we let S r (x) = 
min(<5(x), r) and then we require that 

f 5ida^r n - 1+t , t>-l + e (2.2) 

J A r (x) 

f <5* dy « r n+ \ t > -2 + e. (2.3) 

In our main results there will be a restriction on e. The condition on e depends 
on the L 9 -index in the reverse Holder inequality in Theorem 13.11 and we refer 
the reader to section [4] for more details. For the moment, we observe that given 
M we will find a positive value of e for which the conditions (|2.2p and (|2.3p 
allow us to solve the mixed problem. While the conditions (12.2p (12.3P may 
seem rather mysterious, they are closely related to the dimension of the set A. 
The work of Taylor, Ott, and Brown |33^ Lemmata 2.4, 2.5] shows that the 
conditions (I2.2p and (|2,3p follow if we assume that the set A is of dimension 
n — 2 + e. In the setting treated in [26J, we have e = 0. As the conditions 
(12. 2p and (]2.3p become more restrictive as e decreases, the results of this paper 
always apply in the domains considered in [26]. 

Example. We close this sub-section by giving an example of a domain that 
satisfies our conditions and illustrates that even in two dimensions, the set D 
can be fairly complicated. Thus let f2 = {x : \x\ < 1} be the disk and let D be 
the set U^ 1 {(cos0,sin6') : 7r/2 2fc+1 < 6 < vr/2 2fc }. It is not difficult to see that 
this domain will satisfy (f2TT]l and (j2~2j) and (f2~3j) with e = 0. 



2.2 Function spaces 

We will consider L p spaces with respect to Lebesgue measure on domains in 
C R n and with respect to surface measure, a, on the boundary of f2, <9f2. 
These will be denoted L P {Q) and L p (dQ), respectively. 

We let W l,p {Q) denote the standard Sobolev space of functions having one 
derivative in L p (Vt) with the norm 

IMIw 1 ^) = (J \Vu\ p + r Q p \u\ p dy 
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at least for 1 < p < oo. The factor of vq guarantees that the two terms in the 
norm have the same homogeneity when rescaling. Our functions will generally 
take values in R n , though we will not indicate the range in our notation. 

12 19 

For D C <90, we define WjJ (Q) as the closure in W ' (O) of the functions 
which are smooth in the closure of Q, and which vanish in a neighborhood of D. 
We let Wp 1,2 (Cl) denote the dual of W^ 2 (Q). As noted all domains in this paper 
will be Lipschitz and thus we have the trace operator which is a continuous map 
from W 1 ' 2 ^) into L 2 (dQ) and extends the operation of restricting a smooth 

1 /2 2 12 

function to the boundary. We let W D ' (dfl) denote the image of w^(n) 

/q 2 1 /2 2 

under the trace map and then W D (dVL) denotes the dual of W% ' (dSl). 

We will also need to consider Sobolev spaces on the boundary. For a function 
u in C°°(Cl) we define a family of tangential derivatives 

du du du 

a — = u ia v iTT~ ' hJ = ^,---,n,t^j 

OTij OXj OXi 

where v denotes the unit outer normal to the boundary. For a function u which 
is smooth in Q, we can define the tangential gradient by 

V t u a = Vu a - vVu a -v, a = l,...,n. 

Note that we have 

IV.7/1 2 » V 

dn. 



l<i<j<n J 



and for 1 < p < oo, the Sobolev space W 1,p (dfl) is the closure of the boundary 
values of smooth functions in the norm 

\H\w^(dn) = (\Ntu\\ p LP{dn) +r~ p \\u\\ P LP(m) ) 1 / p . 

Before defining Hardy spaces on the boundary, we recall the definition of 
the spaces of Holder continuous functions. For K a compact subset of Cl and 
< a < 1, we define the Holder space C a (K) to be the collection of functions 
/ on K for which the norm below is finite 

Wf\\c(K) = ™p '^"jy 1 + r " Q sup |/(x)|. 

x+y \X-y\ K 

If Q. is a Lipschitz domain in dimension n and N C d£l is a closed set, 
we introduce the Hardy space H P (N) for (n — l)/n < p < 1. We say that 
a is an atom for dQ if a is supported in a boundary ball A r (x), a satisfies 
ll a ll-L°°(dn) ^ l/o"(^r(a;)) and J A , . ado = 0. We say that a is an atom for N 
if a = a | at for some a, an atom for dO,. As we shall see these atoms are building 
blocks for the spaces H P (N) which arise naturally in the study of L p -mixed 
problem for p < 1. We note that our atoms are normalized for the space H 1 . 
When we introduce H p spaces with p < 1, we will need to introduce powers 
of <r(A r (x)) into the norm to compensate for the fact that the atoms are not 
normalized in H p . 
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For p satisfying (n — l)/n < p < 1, we say that / lies in the Hardy space 
H P (N) if there exists a sequence of atoms for N, {aj}JL 1 and coefficients 

{Xj}f =1 C (0,oo) so that / = ET=i X J a J and ET=i XP j (T ( A rA x j)) 1 ~ P < 00 
where the atom aj is supported in a boundary ball A r (xj). We define a quasi- 
norm for p < 1 by 

H/IU,(iV)=rnf ^A^(A r3 (x J )) 1 

where the infimum is taken over all representations of / as a sum of atoms. 
When p = 1, the sum converges in L 1 and we have H 1 ^) C ^(N). However, 
for p < 1 elements of H V {N) may not be functions. Rather they are defined 
as linear functionals on a space of smooth functions. It is well-known that 
elements of H p (d£l) give continuous linear functionals on the space C a (d£l) 
where a = (n — l)(l/p — 1), for 1 > p > (n — l)/n. It is a straightforward 
consequence of our definition that elements of H P (N) give continuous linear 
functionals on C£>(<90), the collection of functions in C a (d£l) which vanish on 
D. We observe that it is an immediate consequence of the definition that an 
element / from H P (N) has an extension to H p (d£l) with the bound ||/|| j^p(asi) ^ 
C||/||# P (jv) for any C > 1. One consequence of our main theorem is that there 
is a bounded linear extension operator from H P (N) to H p {d£l). However, it 
would be interesting to develop a better understanding of these spaces. 

We define Hardy-Sobolev spaces H 1,P {D) as follows. We say that A is a 
1-atom for dQ, if A is supported in a boundary ball A r (x) and || V^l^oo^) < 
l/a(A r (x)). We say that A is an atom for D, if A is the restriction to D 
of an atom for d£l. Finally, if p satisfies (n — l)/n < p < 1, we say u is in 
H 1,P (D) if there is a sequence of atoms {Aj} for D with each Aj supported in 
a boundary ball Aj = A r \xj) and a sequence of non- negative real numbers 
{Xj}JL 1 with J2j Xjcr(Aj) 1 ^ p < 00 and u = J2j X jAj. We define a quasi- 
norm on H^ p (dn) as \\u\\ H i, P{dQ) = mf(J2j X p j a{A j ) 1 ' p ) 1 / p where the infinum 
is taken over all possible representations of u. It is well-known that one may 
define atoms using an L l space 1 < t < 00 instead of L°°. Thus, if we fix 
t with 1 < t < 00 and replace the condition ||a||L°o(A r (2)) < c(A r (x)) _1 by 
||a||x,t(A r ( x )) ^ cr(A r (x))^ 1 ^ t , then the resulting Hardy space is independent of 
t. A similar result holds for the Hardy-Sobolev space. The work of Mitrea and 
Wright [241 sections 2.2-2.3] provides an exposition of these facts in the setting 
of the boundary of a Lipschitz domain. 
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2.3 The weak formulation of the mixed problem 

Next we recall the Lame operator L and give several well-known estimates for 
solutions of the Lame system, Lu = 0. The operator L will be written as 

d d 



s 



Here and throughout this paper, we will use the summation convention and 
thus we sum on the repeated indices i, j and /3 in the above expression. We 
will consider the traction operator at the boundary 

du a _ ij dv^ 
dp 1 01/3 dxj 

The coefficients a^a we consider are given by 

a «/3 = a a/?( s ) = f^ijfiap + sSipdja + (A + p - s)S ia 6jp (2.4) 

where 5{j is the Kronecker symbol, A and p are the Lame parameters and s is 
an artificial parameter that has been introduced in several different settings for 
mathematical reasons. Note that 

d d 

Thus, the differential operator L does not depend on the value of s. However, 
the traction operator du/dp does depend on s. 

We will require that our operator satisfy the ellipticity condition 

a%e a e p >M-^-±f\\ eeR nxn . (2.5) 

This will hold if s, A, and p satisfy 

s € [0, p], p> 0, and 2p + nX > 0. (2.6) 

We note several consequences of this ellipticity assumption (|2.6p which come 
under the name of Korn inequalities. First we have the global coercivity esti- 
mate 

c / o iv.r*</ o «&*£*£*, uew?m m 

This will hold if D is a non-empty open subset of the boundary. The constant 
will depend on Q, and D and the constant will be bounded away from zero for 
a family of operators which satisfy (|2,5p for a fixed M. 

For u G W 1,2 (£l), any local domain Cl p (x), and any constant vector c G R n , 
we have 

/ \Vu\ 2 dy<C(M) [ a % d ^-^+p- 2 \u-c\ 2 dy. (2.8) 
Jn p (x) Ju p (x) p d yj dy % 

The constant depends only on M and the estimate does not depend on the 
properties of D and, in fact, we do not require that u vanish on D for the 
estimate (|2.8j) to hold. A proof of (|2.8p may be found in the monograph of 
Oleinik, Shamaev, and Yosifian |25l Theorem 2.10]. Throughout this paper, 
we will consider the Lame operator (and its generalizations) with coefficients 
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given by (|2.4p and satisfying (|2,6p . As a consequence of these assumptions, we 
obtain the estimates (i23]l . (f277|) . and (j2~H . 

We are ready to give the definition of a weak solution of the mixed problem. 
Given / G W^ 1,2 (Q) and /at G VF~ 1 / 2 ' 2 ((9il), we consider the boundary value 
problem 

Lu = /, in 

u = 0, on D (2. 



§ = /at, on iV. 



We say that u is a weafc solution of \2. 9\) if u G W^ 2 {VL) and 



n 



a as^— d V = (fx, <f>)an - {f, 0)n, </> G W^' (0). 



We are using (•, -)q and (•, -)qq to denote the pairings of duality on W D 1,2 {£1) x 

2 1/2 2 1/2 2 

W^' (0) and W D 1 * (dU) x W£ ' (dQ), respectively. As we assume that the 
coefficients of L satisfy (|2.5p and we assume D is a nonempty open set of 
d£l, we have the coercivity estimate (|2.7p and a Poincare inequality (|2.10p . If 
u G W^'^fi), then we have that 

IMIiP(n) < Cr \\Vu\\ LPin y (2.10) 

Thus the existence and uniqueness of solutions to (j2.9j) is a straightforward 
consequence of standard Hilbert space theory. The solution of (12. 9p will satisfy 

W Vu \\w^ 2 (n) ^ c (\\f\\w-^(n) + ll/iv|| w -iA2 (an) ). 

Now that we have introduced the traction operator, d/dp, we describe the 
sense in which we take boundary values of du/dp when p < 1. As the boundary 
values may not be a function, it is not enough to ask for non-tangential limits 
a.e. Thus for p < 1, and / G H p (dQ), we say that du/dp = / on dtt, if we have 

lim / ^{^fda={f^) d n, ^C Q p (2-11) 

where £lk C SI is a family of approximating domains as constructed by Verchota. 
If /tv lies in H P (N), we say that du/dp = /at on TV" if we have (|2.1ip for which 
lie in C%{d£l). 

One technical point about our argument for the mixed problem is that 
we do not have solutions in smooth domains as a starting point. A common 
strategy in studying the Dirichlet and Neumann problems is to approximate a 
Lipschitz domain by a sequence of smooth domains, prove a uniform estimate 
in the approximating domains, and take a limit. It is an interesting question to 
find an approximation of the mixed problems by problems which have smooth 
solutions. 

Finally, we make a note about the constants in our estimates. Many of 
our estimates are of a local nature and hold on scales r with < r < ro and 
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with a constant that depends only on the parameter M that appears in the 
definition of a Lipschitz domain, the corkscrew condition, the constants in the 
estimates (|2.2p . and (|2.3p . and any L p indices that appear in the estimate. We 
will say a constant depends on the global character of il if the constant also 
depends on the collection of coordinate cylinders appearing in the definition of 
the Lipschitz domain and the constant in the Korn inequality (|2.7p , 

3 A reverse Holder inequality for weak so- 
lutions 

The first step of our argument is to show that if u is a solution of (|2.9p with nice 
data, then Vit in L g (£l) for some q > 2. We use the reverse Holder argument of 
Gehring [14J and Giaquinta and Modica |16j . The argument is similar to one 
given in our joint work with Taylor [261 E2], but a few changes are needed since 
the operators involved may not be strongly elliptic. The next Theorem gives 
a precise formulation of our result. We note that the proof of this result only 
requires the condition (12.ip on D. 

Theorem 3.1 Let Q C R n be a Lipschitz domain and suppose that D satisfies 
| QOP and that L is the Lame system with coefficients satisfying \2. 6\) . Suppose 
that u is a solution of \2. 9\) with f = and fx a function. We may find qo > 2 
so that if 2 < q < qo and fx & L q ^ n ~ 1 ^ n (N), then for any local domain Q r (x) 
we have 




(3.2) 

The value of qo depends only on M. The constant in this estimate depends only 
on M, the dimension n, and q. 

In the above theorem and throughout this paper, we use f- E f to denote the 
average of / over a set E and we adopt the convention that the average over 
the empty set is zero. 

We introduce the notation u x ^ r which will be useful for the arguments of 
this section. For u a function defined on $7, x E ft and r S (0, lOOro), we let 
u x ,r = if 5(x) < r and u x , r = ,su(y)dy if 5(x) > r. We observe that 
one useful feature of this definition is that if r\ is in C°°(Q) and 77 = outside 
Q r (x), and u S Wjj (fl), then rj(u — u x>r ) lies in W^ 2 (Cl). 

We give several formulations of the Sobolev inequalities that we will use in 
the sequel. The proofs of the estimates (|3.3p and (|3.4p are given in Appendix [Al 
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These take advantage of the corkscrew condition for D (|2.1j) and the definition 
of u xr . Let p and q satisfy 1 < p < n and \jq = 1/p — 1/n. Then for x € 
and < p < r < lOOo, we have 

For 1 < p < n, we define q by 1/q = n/(p(n — 1)) — l/(n — 1). We observe that 
for all local domains Q p {x) and Q, r {x) with < p < r, we have 

(7 l»-^*Y Vf »?£'' (/ |V<*) 1/P . (3.4, 

We turn to the proof of the reverse Holder inequality in Theorem 13.11 The 
proof will require the following auxiliary function 

P r f(x)= sup — j- / \f\da. 

0<s<r s JdnndQ a (x) 

When we apply P r to the traction data /at we will assume that /at is extended 
to be a function on dil by setting /jv(x) = on D. As observed in [26], for 
1 < p < oo, q = pn/(n — 1), we have 

WPrfh^n) < C\\f\\ LP{dn) . (3.5) 

Proof of Theorem \3.1l We let u be a weak solution of (|2.9p with / = and /at 
a function. We fix a; € d£l, p and r with r /2 < p < r < r^ and set e = (r — p)/2. 
We claim that with p satisfying p > 2(n — l)/n, we have 




(3.6) 



Observe that the last term is bounded by P2 r (\fj^\ p ) l ' p (x). Theorem 13.11 will 
follow from (|3.6p using the argument in Giaquinta [15, p. 122], and the estimate 
([33]). 

Thus, we turn to the proof of (|3.6|) . We let rj be a cutoff function which is 
one on £l p (x) and zero outside Q p+e (x) and observe that v = if(u — u XiP+e ) G 
Wjj (Q). We let E = J n + ^ \u\\Vu\rj\Vr]\ dy. Using the product rule, the 
local ellipticity assumption (I2.8p . and that u is a solution, we obtain 

f V 2 \Vu\ 2 dy = [ \V V (u-u x , p+e )\ 2 dy + CE 

JQ. p+t (x) JQ P+I :(x) 



f&vP da + CE. 



an 
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We apply Holder's inequality, Young's inequality with e's, and the boundary 
Poincare inequality (|3.4p and obtain for any 7 € R that 



/ f N v^da\<r^([ \f N \ p da] 
Jan \Janndn p+e (x) J 



\ 2/P 

(3.8) 



CV~ 7 f f \ 2 /* 

+ (W^U^. |v " l ' <!y j ' 

where p > 1 and p and t are related by 1/t = l/n + (n — 1) / (np 1 ). We return to 
(|3.7p and use (|3.8j) to estimate the boundary term. To handle the error term 
E, we use Cauchy's inequality with e's to subtract the term Vu in E from the 
left-hand side of (|3.7p and conclude that 

f rf\Vu\ 2 dy<C(r- 2 [ \u - u x>p+t \ 2 dy + r 7 [ \f N \ p da) 
Jn p+C (x) Jn p+C (x) \Jdnndn p+i (x) J 



r 



-7 



v 2/t 

Vu|*dy ) 



(1 _ p/r) 2n-2 \J Up+2e{x) ' ' J 



(3.9) 



We use (|3.3p to estimate the first term in (j3.9j) . choose 7 so that n + 7 = 2n/p 
and hence n — 7 = 2(n — 1) /p, divide by r n , and recall that 2e = r — p to obtain 

/ \ 2/p 

/ r/ 2 |Vu| 2 dy<C[ l-^j / \f N \ p da) 

Jn p+i (x) \ r Jdn r (x)dn J 



){n+2)/n / \ 2/t 

+ (/ |Vn|*dy )] 
y7n r (a:) y 



+ (l-(p/r))2-2 ( ^ rW 

We let s = max(t, 2n/(n+2)) and observe that our condition that p > 2(n— l)/n 
implies that t < 2. We choose 9 so that 1/s = (1 — + 0/2 and use Holder's 
inequality and Young's inequality to obtain that 



\ 1/2 

|Vu| 2 dy] < C((l- p/r)- 2 ^ -f |V«|d 2 / + P r (|/n 1 / p (a;)) 



+l(~f \Vu\ 2 dy) 
IU r (x) J 



2 

Now a standard argument as in Giaquinta [TU p. 161] gives f)3.6f) . I 

4 Atomic Estimates 

Our next step is to estimate the gradient of solutions at the boundary and 
eventually to obtain estimates for the non-tangential maximal function (Vit)*. 
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It is this part of the argument that requires the integrability conditions on 5 
given in (J22J and (l2~3j) . 

Our main result is Theorem 14.11 A key step in the proof, Lemma 14.101 uses 
estimates for the Green function from the work of Taylor, Kim, and Brown [32]. 
We are only able to prove these estimates in two dimensions. 

Throughout this section, we assume that ft is a Lipschitz domain in R 2 , 
that D C dSl is an open set satisfying (|2.ip . and that L is the Lame operator 
with coefficients satisfying (|2.6p . If qo is as in Theorem 13. 11 we choose e so that 
1 < qo(l — e) / (2 — e) and require that Q and D satisfy (|2.2p and (12. 3D for this e. 

Given a boundary ball A r (x), we let So = S r (x) = B r (x)ndQ and for k > 1, 
we let Efc = 5 , 2 fe r ,(x)\S , 2 fc-i r (x). We will use this notation in the following proof 
and in a similar argument in Appendix [Bj 

Theorem 4.1 If u is a weak solution of the mixed problem \2. 9\) with f = 
and traction data /at = a an atom supported in A r (x), then for p with p < 
qo(l — e)/(2 — e), we have 

ll(Vn)*|| LP(Efe) < C2- k ~<(2 k r)^-\ k > 0. (4.2) 

Here, 7 is the exponent of Holder continuity for the Green function for our 
mixed problem (see f3Ej or the estimate \4. 13\ ) below). In addition, Vn has 
non-tangential limits a.e. on dfl. 

Let pi = 1/(1 + 7) then for p\ < p < 1, we have 

\\(Vuy\\ LP{dn) < Ca(A r (x)f P -\ (4.3) 

The constants in these estimates depend on p, M and the global properties 
of SI. 

We recall the following results of Dahlberg, Kenig, and Verchota [TTJ which 
are a consequence of the existence theory for the L 2 -traction and regularity 
problems. These results do not depend on D and hold in all dimensions. To 
state these estimates, we will use the truncated non-tangential maximal function 
defined by 

v*(x) = sup |Vf(y)|, x £ dQ. 

yeT(x)nB r (x) 

Lemma 4.4 ( [11] ) Let x G dQ and r S (0, 25ro). Suppose that u is a solution 
of Lu = in Q^ r (x) and u lies in the Sobolev space W 1 ' 2 (Q l i r (x)). 
Ifdu/dp e L 2 (A Ar {x)), then 



[ (Vu)fda<C 

JA r (x) 

IfVtu lies in L 2 (A Ar (x)) , then 





du 




I 


dp 


da + - f 


\ iA 4r (i) 


r Jn 4r (x) 



|Vii| 2 dy 



f {Vu)f da < C I [ \V t u\ 2 da + - f \Vu\ 2 dy 
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In each case, Vu has non-tangential limits a.e. on A r (x). The constants 
in these estimates depend only on the dimension and M . 



These estimates may be proven using arguments from [TT]. See Mayboroda 
and Mitrea |21| for results in two dimensions. 

We give a notion of Whitney decomposition on the boundary of a Lipschitz 
domain. To pass from R n to a more general set, we sacrifice the property 
that an open set is decomposed into disjoint sets and only require that each 
point lie in at most finitely many elements of the decomposition. This seems 
to be enough for our applications. Given an open set O C dQ, a Whitney 
decomposition of O is a collection of boundary balls {Aj = A Tj (xj)}JL l so 
that 1) O = UjAj, 2) for a constant C\ that may be chosen as large as we 
like, we have C X Tj < dist(A i ,<90 \ O) < 2C\Tj, 3) \o < J2jXA, < c n xo, 4) 
Ylj Xn 4r ,(xj) < C- We leave it as an exercise to construct this decomposition. 

Lemma 4.5 Suppose u is a weak solution of the mixed problem \2. 9\) with /at 
in L 2 (N) and / = 0. Then for x 6 dQ, r £ (0, 50ro) and p € R, we have 

[ (Vuy cS *8 l -Pdo-<c[ f \f N \ 2 5 1 r ~ p d*+ [ \Vu\%»dy). 

JA r (x) \JA 2 r{x) JQ 2r {x) J 

The constant in this estimate depends only on M, the dimension n, and p. 

Proof. We construct a Whitney decomposition of <9il \ A and let {Aj = 
A rj (xj)} denote the boundary balls from the Whitney decomposition which 
intersect A r {x). This gives a family of boundary balls so that for each j, 
A r (x) \ A C UjAj C A 2r (x), A 2rj (xj) C N, or A 2rj (x i ) C D, rj w S r (y) for 
y € fl^rj {x), and ^2xn 4r .(x) — C(M,n). We apply the estimate of Lemma 14. 4^ 
sum on j, and use the properties of the Whitney decomposition to obtain the 
Lemma. I 

The next lemma takes us from estimates in a weighted L 2 space to estimates 
in an unweighted LP space. Recall that e appears in our hypotheses (|2.2p and 
(|2.3j) and is required to satisfy 1 < qo(l — e)/(2 — e) with qo as in Theorem 13.11 



Lemma 4.6 Let u be a weak solution of \2. 9\) with /at in L°°(N) and / = 0. 

Let pq = qo(l — e) / (2 — e) and suppose that 1 < p < pq. For each boundary 
ball A r (x) with < r < 50ro- 

Vp 

<C{4 \Vu\dy + ||/jv||L°°(A 2r (z)))- 

JQ.2r{x) 

The constant in this local estimate depends on M and the constants in K2. 2\) 
and (KM>- 
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Proof. We choose p and q with 1 < p < q(l — e)/(2 — e) < qo(l — e)/(2 — e). 
We then fix p so that (2 - e)(l - (2/q)) > p > (2 - e) - 2(1 - e)/p. We apply 
Holder's inequality with exponents 2/p and 2/(2 — p), the property f)2.2|) . the 
lower bound for p and Lemma 14.51 to obtain 




(4.7) 

To estimate the integral over £l2r(x), we use Holder's inequality with expo- 
nents 2/q and qj (q—2), the property (|2.3p . the upper bound for p, and Theorem 



13. II to obtain 




q(n-l) 



(4.8) 

The Lemma follows from (|4,7p , (|4.8p , and a simple covering argument to replace 
4r by 2r on the right hand side of the estimate. I 



Before we proceed to the proof of Theorem 14.11 we recall a version of the 
energy estimate for solutions of the weak mixed problem, f)2.9|) . 

Lemma 4.9 Let u be a solution of the weak mixed problem \2. 9\) with data 
/at = a, an atom for N and f = 0. The solution u satisfies 

[ \Vu\ 2 dy < C. 
Jn 

The constant C depends on the global character of CI. 

Proof. According to ([33]), we have that W^ /2 ' 2 (dn) C BMO(dQ), the space 
of functions of bounded mean oscillation on dd. Thus if a is an atom for TV the 
map u — > f dn a^u 13 da lies in Wp 1,2 (Q). The estimate of the Lemma follows. I 
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We are only able to prove the following Lemma in two dimensions. Finding 
an appropriate substitute for this Lemma is the main obstacle to studying the 
mixed problem for the Lame system in higher dimensions. 

Lemma 4.10 Let a be an atom for N that is supported in a boundary ball 
A r (x). We may find an exponent 7 > which depends only on M so that if 
u is a solution of the weak mixed problem \2. 9\l with fjy = a and f = 0, then 
with p as in Theorem \4-l\ we have a constant C so that 



The constant in the estimate depends on the global character ofQ. 

This Lemma gives the estimate of Theorem 14.11 for Vu. Additional work is 
needed to obtain the estimate for the non-tangential maximal function. 

Proof. We fix a boundary ball A r (x) and let a be an atom for iV that is 
supported in A r (x). We let u be a solution of the weak mixed problem (12.9P 
with /tv = a and / = 0. We apply Lemma I4T61 on a boundary ball Ajj(y) and 
use the normalization of the atom to obtain that 



When n = 2, the Cauchy-Schwarz inequality and the energy estimate of Lemma 
14.91 gives that for any local domain £ln(y) 



and if we set y = x, we obtain the estimate (|4.11j) but with a constant that 
depends on k. 

We fix ko so that we have dist(A r (x), £fc ) ?» r and establish (|4.1ip for 
k > k$ with constant that does not depend on k. Towards this end, we consider 
a boundary ball A#(y) with Ag#(y) n A r (x) = 0. Since a = on A2_r(?/), the 
estimate of Lemma 14.61 implies that 




k > 0. 



(4.11) 





Together, these observations imply that 





The local ellipticity condition (|2.8p implies a Caccioppoli inequality 




(4.12) 
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provided that u is a solution of (|2.9p with /jv = on A4#(y) and / = on 
^4fl(y)- Next we note that the Green function estimates of Taylor, Kim, and 
Brown [32] and the arguments in [26] imply that we have a constant C and 
exponent 7 > so that 

Hv)\ <c(-^— V, i/en\fl 2r (x). (4.13) 

V 1^ — 2/1/ 

Combining the estimate of Lemma 14.61 the Holder inequality, (|4.12j) . and (|4.13p 
gives that 

C ( r \ 7 



\Vu\ p da\ < — 



We may cover by boundary balls with radius comparable to 2 k r and use the 
above estimate to obtain the estimate of the Lemma for k > ko- I 

Now we are ready to give the proof of Theorem 14.11 



Proof of Theorem \4.1\ We let <fi E C°°(ti) be a smooth function which is zero 
outside ^4 r (x) and is zero in a neighborhood of A. We let T a P be the matrix 
fundamental solution of the Lame operator in R 2 , thus (LP 7 )" = <5q. 7 <5o where 
<5 a7 is the Kronecker delta and 5q is the Dirac delta measure at 0. We claim 
the following representation formula for derivatives of u, with a,k = 1,2, 



du a f d ■■ dT^ du a 

m ^ x)= L^ l ^ {x -- )ck w k dy 



dT^ , , du a dT^< , x , du a 



x 



dn ^ dyj dy k ali dyj % 

ij dT^ d(p du a ij dT^ d<p du a 

(4.14) 

Of course, the integral to the right of the first equals sign should be interpreted 
as the Dirac delta acting on a smooth function. If u is smooth, the proof of 
(|4.14p is accomplished by integration by parts. We begin by moving the d/dyi 
derivative, then the d/dy^ derivative, and finally the 8/ dyj derivative. Since 
Lemma [4.41 gives non-tangential maximal function estimates for u and Vu away 
from A and we assume that <f> is supported in a coordinate cylinder and is zero 
in a neighborhood of A, standard limiting arguments allow us to obtain the 
representation formula (|4.14p without the assumption that u is smooth in f2. 

Next, we remove the restriction that <f> vanish near A. Let ip T be a function 
which is one when 5(x) > 2r, zero if 5(x) < r, and satisfies \d a vp T /dx a \ < 
C a /rl a L Replace <fi in (|4.14p by t[j T <fi where <fi is a smooth function which is zero 
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outside Q,4 r (x), but is not required to vanish on A. Using Holder's inequality, 
the terms involving derivatives of ip T in (14.14p have the upper bound 



f \VT{x--)\\Vu\\<j)Vip T \dy< — ( I \Vu\ q dy] |A T , r | W (4.15) 

JA Ti r r WA Tir I 

where A T)r = Q/i r (x)n{y : r < 5(y) < 2r} and q is chosen to satisfy (2 — e)/(l — 
e) < q < qo. From the estimate (|2.3p . we obtain the upper bound 



|A T , r | < Cr a / 5- a dy<C r T a 

provided < a < 2 — e. Our choice of g implies </ < 2 — e so we set a = q' . 
Then the estimate (j3.2j) for Vti implies that the right-hand side of ()4.15j) tends 
to zero with r. From the estimate of Lemma 14.10} we have Vu is in L p (d£l) 
for some p > 1, thus we may use the dominated convergence theorem and let r 
tend to zero in the boundary terms of (|4.14p . This removes the restriction that 
(f> vanish in a neighborhood of A. 

To establish ()4.2|) . we begin with the representation formula ()4.14p . We let 
<p run over a partition of unity and sum to obtain the representation 

du a , , f a dT^ . . , du a du a ^ dT^ . . H du a , , 1 n 

The divergence theorem implies that 

^ a a«l^ d(J = ' a = 1,2 (4.16) 

'on 

da = 0, i,k,P = 1,2 (4.17) 



where the techniques used to prove (|4.14p help to justify the application of the 
divergence theorem. We may use the theorem of Coifman, Mcintosh, and Meyer 
[6] , the observations (j4.16M.17]) , the size estimates of Lemma I4.10j and the 
notion of molecule (see [7]) to establish the estimates (|4.2p . For p > 1/(1 + 7), 
the estimate (|4.3p follows easily from (14. 2p . I 



5 Uniqueness 

We take advantage of the restriction to n = 2 to give a different proof of 
uniqueness than found in [261 [33] . The argument uses a version of the Hardy- 
Littlewood-Sobolev theorem on fractional integration. The proof is adapted 
from one found in [3] and [24\ Lemma 11.9] where a result similar to Lemma 
15. H is established in dimensions n > 3. The result for n = 2 is more technically 
involved. 
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Lemma 5.1 Let L be the Lame operator with coefficients satisfying \2. 6\) . Let 
O C R 2 be a Lipschitz domain and suppose that u is a solution of Lu = with 
(Vu)* G L p (<90). If p < 1, q is defined by 1/q = 1/p — 1, and K is a compact 
subset of O, then we have 

IKIIw(dfi) < C(\\(Vu)*\\ LP{m) +ro /p ~ 1 sup|u|). 

The constant in this estimate depends on p, K and the global character o/O. 

For this argument it will be convenient to use S r (x) = B r (x)DdO, for x € <90 
in place of A r (x) and to define the Hardy-Littlewood maximal operator on the 
boundary by 

Mf{x)= sup 4 |/| da. 

r>0 Js r (x) 

Proof. We let {0&} be a sequence of approximating domains with 0& C O as 
in Verchota [36] . We will show that we have 

ll«*IU«(«i fc ) < c (\\(^ u )*hp(dn k ) +rl /q sup\u\) (5.2) 

K 

where p and q are as in the statement of the Lemma and the constant C is 
independent of k. Letting k tend to oo gives the result for 0. The advantage of 
working on the subdomain 0& is that we may assume that 1 1 1 1 (<9fi fc ) ^ s finite. 
In what follows, we drop the subscript k and assume that ||u*||L9(9f2) i s finite. 

To begin we observe that if y £ T(x), for some x € <90, then for e > 0, and 
/3 a multi-index, we have 




Here, S(x,y) = {z : \z — x\ < (2 + 2a)d(y)} n <90 with a as in the definition 
of the approach regions T(x) and recall that d(y) = dist(y, 90). To establish 
(|5.3p . choose x € 90 so that \x — y\ = d(y). We claim that if \w — y\ < 
min{a/4, ct/(2 + 2a)}d(y) and \z — x\ < aci(y)/4, then u; € T(z). To establish 
this claim, first note that d(y) < d(w) + \w — y\ < d(w) + ad{y)/(2 + 2a) 
and hence d(y) < d(w)(2 + 2a)/ (2 + a). Using the triangle inequality and 
our choices for z , x , and u>, we have |z — io| < [io — y| + |y — x\ + |x — z\ < 
(1 + a/2)d(y) < (1 + a)d(w). Also, if z € S a d(y)/i(x), then by the triangle 
inequality |z — x| < \z — x\ + \x — y\ + y — x\ < (2 + 2a)d(y) and we conclude 
that S ad ( y y 4 (x) C S , ( 2 +2a)d(j/)( :c )- Tnus > with r = min{a/4, a/(2 + 2a)}d(y), 
we may use interior estimates for derivatives of u to obtain 



(y)| < „ Mo, sup |u(w)| < - , , 4 u* e da 



w£B r (y) JS(x,y) 
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The second inequality follows since B r {y) C T(z) for all z in S a M y \u{x) and 
hence S a M y \{x) C S(x,y). We have established (I5.3p . 

To begin the proof of (15. 2j) . we let {Zj = Z^a^)}^ be the covering of 
90 by coordinate cylinders guaranteed by the definition of a Lipschitz domain 
and recall that we are in R 2 . We let Tj = {(2/1,2/2) : \yi — < ^o>y2 = 
£j,2 + (4M + 2)ro} denote the "top" of each cylinder Zj and then let if' denote 
the compact set UjTj. Applying the estimate (|5.3j) to with /3 = and e = p 
gives that for any pair K and J of compact subsets of 0, we have 

sup|u| < sup|u| + Crl~ 1/p \\(Vu)*\\ LP(m) . (5.4) 
J K 

The heart of the proof is to establish for e > 0, 8 G (0, 1), and x G 90, that 

<C(M(u* e )(x)^/ ((Vu)*)(x) e + sup|u|+r o 1_ ^KV«ril^(an)) (5-5) 

where Ig is the fractional integral defined by 

Ief(x)= [ f{y)\x-y\ e - l da. 
Jan 



We accept the claim (|5.5|) for the moment and complete the proof of the 
Theorem. We raise both sides of (I5.5|) to the power q, integrate over 90, and 
apply the Holder inequality with exponents 1/6 and 1/(1 — 6) to obtain 

IKIIwOO) < C'(||A<(w«)||£7^ ) ||^((V«r)||£» (OT) H-rJ / «8up |tt|+|| (Vtt)* 

We choose e so that q/e > 1 and hence the maximal operator A4 is bounded 
on L«/ e (90), choose with < < p so that /<, : LP/ e (90) -> L^ e {dn) with 
6/q = 6(l/p — 1). With these choices, we may use Young's inequality, our a 
priori assumption that u* G L 9 (90) and (j5.4j) to obtain the Theorem. 

Finally, we establish the key estimate f)5.5f) . We fix x G 90 and consider 
y G T(x) and suppose that ?/ lies in a coordinate cylinder Zj. If a is large 
enough, we have that the line segment {y + te 2 ■ y 2 < t < ro} lies in T(x). 
(As usual, we are using the coordinate system for the coordinate cylinder Zj to 
state this condition.) Thus, we may use the fundamental theorem of calculus 
to write 

f r ° du 

\u(y)\< \-— (y 1, t)\dt + sup \u\. (5.6) 
J y2 oy 2 k> 

We apply the estimate ()5.3[) to obtain the bound 

dy 2 d{y 
Also, if we apply estimate (|5.3|) to Vu, then we obtain 



(yi,*)! 1 -* < ^ ^MA, 



\Vu{ yi ,t)\ e <C J {Vu)* e da. 

Js(x,(yi,t)) 
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We write \du/dy2\ 1 9+9 in ()5.6[) and use the two previous estimates to obtain 

ro 

\^{Vx,t)\dt<CM{u*'){x)W/< 
y2 dy 2 

x I ° ,„ 1 ., 2 g / (Vu)* e dadt ( 5J ) 

< CM(u* e )(x) {1 ~ 9y€ I e ((Vuy 9 )(x). 

The last step uses Tonelli's theorem to change the order of integration. The 
estimates (|5.6p and (|5.7j) give a bound for the supremum of u in n (UZj). 
Estimate (|5.4p allows us to estimate the supremum of u in \ UjZj. Together 
these two observations give (15, 5p . This completes the proof of the estimate 
(|5.5p and hence the Lemma. I 

We are ready to give uniqueness for the L p -mixed problem. 

Theorem 5.8 Suppose f2 C R 2 is Lipschitz domain, f2 and D satisfy \2. 
pLi| ), and \2. 3\) with ^q(1 — e)/(2 — e) > 1, anc? L zs the Lame system with 
coefficients satisfying \2. 6\) . There exist p\ < 1 so that if p > p\ and u is a 
solution of M-l\) with /at = and frj = 0, then u = 0. We assume that the 
traction operator exists in the sense of non-tangential limits if p > 1 and in the 
sense of 112.11]) if p < 1. 

Proof. It suffices to prove the Theorem for p < 1. We let a be a solution of 
the mixed problem (jl.ip with /at = and /d = 0. Our argument will proceed 
by duality and requires the existence of a solution when the traction data is 
an atom for TV as given in Theorem 14.11 We let a be an atom for iV and claim 
that 

aPu^ da = 0. (5.9) 



N 

We accept the claim and complete the proof. We first observe that (Vu)* € 
L p (dQ.) and thus Lemma [5.11 implies that u* € L q (d£l) with l/q = l/p — 1. 
Also, it is easy to see that if (Vu)* G L p {d£l) then u has non-tangential limits 
a.e. The claim implies that these limits are zero. From the work of jllj, we have 
uniqueness for the L^-Dirichlet problem for q > 2 and any Lipschitz domain 
and thus we have u = if q > 2. Since we are in two dimensions, we have q > 2 
if p > 2/3, thus we will require that p\ > 2/3. 

To establish the claim (|5,9p . we let v be the solution of (jl.ip with data 
/at = a and fjj = 0. From Theorem 14.11 we have (Vi>)* € L {dQ) for t < p$. 
Lemma O implies that v* G L*(0fi) for 1/g = 1/i-l. Thus if 1/po + l/Pi < 2, 
then we may find i so that 1/t + l/p < 1 and use the dominated convergence 
theorem to obtain 



l im f vP&fd* 

k^oo J m dp 



u p aP da. 

N 
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The estimates for the Green function in section 4 of [32] imply that there is an 
exponent 7 so that the solution v lies in C^(dQ) for a < 7. If a > 1/pi — 1, 
then our assumption that du/dp = in H P (N) implies that 



l im f v PAP d a = 0. 



k- 

While the Green identity gives for every k that 



/ v^-u^fda = 0. 
Jdn k dp dp 



The last three displayed equations imply the claim (15. 9p , 

To summarize the conditions on p±, we need p\ > 2/3, p\ > po/(2po — 1) 
and pi > 1/(1 + 7) to establish uniqueness. I 



6 Existence of solutions 

In this section, we give the details needed to establish the existence of solutions 
for the L p -mixed problem for p in an interval containing 1. When p < 1, we 
take our data from Hardy spaces and for p > 1, the data is taken from LP 
spaces. Given the estimates of Theorem 14.11 the argument is not so different 
from results in previous work of the authors [3] and [26J. 

Theorem 6.1 Let £1 C R 2 be a Lipschitz domain, suppose that D satisfies 
12. \2. 2\) . and \2. 3\) . and that L is the Lame operator with coefficients satis- 
fying the ellipticity condition \2. 6\) . 

We may find an exponent p\ so that if 1 > p > p\, the L p -mixed problem 
has a solution. This means that if /n £ H P (N) and fr> G H l,p (d£l), then we 
have a solution to the L p -mixed problem ([Qp which satisfies 

du 

IMIffL^an) + ll^-lliff(an) + ||(Vn)*||xj>(an) < C(||/2v||hp(jv) + 1 1 Id \\ h 1 * (an))- 

The normal derivative of u exists in the sense of A2.ll]) , u has non-tangential 
limits a.e. on the boundary, and these limits vanish on D. 

Proof. By the results of Appendix |Bl we may assume that fjj = 0. Suppose 
that /at lies in H P (N) and choose a representation of fjy = ^ . Xjaj with the 
atom dj supported in A rj (xj) and Ylj ^o-(A rj (xj)) l ~ p < 2||/jv||hi>(jv). We let 
Vj be the solution of (12. 9p with the traction data /at an atom a,j and / = and 
set um = Z)j=i ^j v j- F r o m Theorem I4.1| we have 

M' 

\\V(u M -u M/ )*\\ p LP(m) <C ^ Aj<r(A P .(a: J -)) 1 - p , M' > M. 

j=M+l 
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From Lemma |5.1| we have that u* M € L q {d£l) for 1/q = 1/p — 1. Thus, the 
estimate of (15.3f) implies that um converges uniformly to a function u on com- 
pact subsets of 0. Furthermore, we have that u is a solution of Lu = in f2, 
||(Vu)*||iP(^) < C||/iv||ffp(jv)) an d Vu has non-tangential limits a.e. on dQ. 
To find du/dp on 90, we first note that as in (I4.16p . we have that for any j 
and k, 



dp 



do = 



and 



/an up Jan k 
The estimates (14. 2D of Theorem 14.11 imply that 



dp 



do = 0. 



,dvj dvj 
~ \Hp(dn) + sup I 1 - - 

k 



q unr(o il} . ii 9 ilHp(an fe ) < C- (6-2) 

To establish (|6.2|) . we observe that the size estimates (|4.2|) and (|4.16|) imply 
that dvj/dp is a molecule as in [7] and thus can be decomposed into atoms. 
Here, we use that the atomic Hardy spaces may be defined with atoms from 
any L space, t > 1, and the space we obtain is independent of t. See Coifman 
and Weiss pi] or Mitrea and Wright [241 sections 2.2-2.3] for more information. 
Thus Y^jLi ^jdvj/dp defines an element of the atomic Hardy space H p (dQ) 
as long as 1/(1 + 7) < P — 1 with 7 as in (14. 2ft . It is also straightforward to 
see that the estimate (|4.2p implies that X/fci defines an element of the 
Hardy-Sobolev space H 1,p (d£l). Now we assume p < 1, fix (f> £ C a (Sl) with 
a = - — 1, and consider the limit 



lim 



an. 



—<t> do- 
op 



OG 



an 



dp 



< lim sup \j / 

fc— >oo ._ 1 J 8Ui 



dp 



< do 



< do 



<%l dp' 



>do 



+ sup 

k 



=A/+1 



an t 



3^ 
<9p 



i do\ + 



an 



0/3 



c/)do\). 



Using the estimate (|4.2[) for the non-tangential maximal function of Vi>j and 
the dominated convergence theorem, we see that the first term on the right of 
this inequality is zero for any M. Since functions in C a , a = 1/p — 1 give rise 
to elements of the dual of H p , the second and third terms are small when M is 
large by our choice of Xj and the estimate for dvj/dp in (|6.2|) . Hence we have 
that du/dp exists in the sense of (|2.1ip and is given by ^j uV j /°P an d we 
have the estimate for du/dp in the estimate of the Theorem. Note that we also 
obtain the estimate when p = 1. Finally, if we restrict 4> to lie in C^(Cl), we 
have that dvj/dp = aj on the set N and thus we have du/dp = /at. I 

Finally, we observe that a real-variable argument of Caffarelli and Peral [5] 
and Shen |28j gives existence for 1 < p < p$. The argument to obtain this 
result is identical to that used in our study of the Laplacian in our earlier work 
with Taylor [26j [33] . 
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Theorem 6.3 Let Q be a Lipschitz domain and let D C <9$7 be a non-empty 
proper open subset of d£l which satisfies \2. \2.2\) . and $2.3\) . Suppose that 
1 < p < po with pq as in Theorem \4-l\ We assume that the coefficients of the 
operator L satisfy $2.6\) . Let f^ G L P (N) and fr> = 0, then we may find a 
solution of the LP -mixed problem 1}) which satisfies 

\\(y u )*\\LP(dn) < C\\fN\\LP(N)- 

The following lemma is a key estimate that is needed to carry out the 
method of Shen [28 . This result may be proved using the techniques in the 
proof of Theorem 14.11 See also the argument in section 6 of our previous work 
[26]. 

Lemma 6.4 Let £1, D, and L be as in Theorem \6.3l Suppose that 1 < p < 
Po = <a>(l-e)/(2-e). 

Lf u is a weak solution of \2. 9\) with data /at in L P (N) and f = 0, we have 
the following local estimate for 1 < p < pq, 

(\ i/p 
/ (Vu); p da) <C{-[ \Vu\dy + U \f N \? da) 1 '?). (6.5) 

JA r (x) J Jn 4r (x) JA 2r (x)nN 

The constant C depends only on M and the exponent p. 



A Sobolev inequalities 

In this appendix, we establish the estimates (]3.3p and (|3.4[) that were used in 
the study of the mixed problem. These results may be found in our earlier work 
[26 1 132 1 [33] . though we do not claim originality. The exposition below serves to 
collect these results in one place and includes an occasional endpoint that was 
missed in our earlier work. 

Let <f) be a Lipschitz function on the unit sphere so that Vt = {y : \y\ < 
r 4>{y/\y\)} is a star-shaped Lipschitz domain of scale r and centered at 0. We 
define a bi-Lipschitz map $ : -Bi(O) — > Q. by 

Hy) = r<f>(y/\y\)y. (A.l) 

Using the argument in |17[ Lemma 7.16] and a change of variables, we may 
find a constant C = C(N,n,p) so that for S C fi, a set of positive Lebesgue 
measure and 1 < p < oo, we have 

f , „ C(M,n,p)r n+1 „ 
\\u - 4 udy\\ LP ( n) < — Wu\\ LP{n) . (A.2) 

J S II 

The only change that is needed from the standard argument is to average u 
with respect to the weight given by the Jacobian of the change of variables. 
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Next we observe that for f2 a star-shaped convex domain with constant M 
and scale r, 1 < p < n, q defined by 1/q = 1/p — 1/n, and S a measurable 
subset of O, we have the Sobolev-Poincare inequality 



uayllz,,^) < ||vw||ii.(n). (A.3) 



To establish fOjl we extend u to R n by a reflection in 90. Thus choose r/ a 
cutoff function which is one when \x\ < (1 + M)r and zero when \x\ > (2 + 2M)r. 
We let 

{ ) \ u(xr 2 <t>(x/\x\) 2 /\x\ 2 )ri(x), x € R n \ £1. 

The Sobolev inequality and the product rule in R n gives 

\\u- 4 udy\\ Lg{n) < \\E(u- 4 udy)\\ Lq{nn) 
J S J s 



< C(p,n)\\VE(u--f «dy)|| L P(R«) 

Js 

< C(l|V«||xp(n) +r _1 ||n- 4 u dy|| LP(n) ). 



Applying the estimate (|A.2j) gives the Sobolev-Poincare inequality (|A.3p . 

The following estimate was proved in [26] under an additional assumption 
that the set D was not too spread out. The proof below is simpler and omits 
that assumption. 

Lemma A. 4 If is a star-shaped Lipschitz domain with scale r and constant 
M , D a measurable subset of dQ, and u S Wp P (Q), then for 1 < p < n and q 
given by 1/q = 1/p — 1/n, we have 

1/9 r n-l / r \ Vp 



The constant C depends on n, p, and M . 

Proof. It suffices to prove this estimate for a function u in C°°(Cl) which 
vanishes in a neighborhood of D. Using the map <E> defined in (jA.ip . we may 
reduce to considering the domain SI = i?i(0). We set A = {y G dBi(0) : <3?(y) € 
D}, then we have 

C{M)a{D) < a(A)r n - x < a{D). 

We will show that if u € C°°(Bi(0)) and u vanishes on A C dB\(0), then for 
1 < p < n 

) 1/p— 1/n / \ 1/p 

<-£n / \Vu\>dy] . (A.5) 
<r{A) \J Bl (o) J 
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As the Lemma follows from the special case where O is the unit ball, i?i(0), we 
only need to prove the claim (1A.5|) . 

We let A = {y € B x (0) :y/\y\ € A} and observe that 



u(y) dy 



< 



n 



A 



\y ■ Vu(y)\dy. 



(A.6) 



To establish ifO]) . we let s be in [0, 1], y 6 ^4 and use the fundamental theorem 
of calculus to write 



1 d / n , 
-u(ty)dt 



y ■ X7u(ty) dt. 



We multiply by s n and integrate on [0, 1] x A to obtain 



u(y) dy 



u(sy)s n 1 dsda 



A Jo 



z ■ Vu(z) dz 



which implies (|A.6p . Given (|AT6|) . the claim §Kj$ follows from (|AT3]> . I 

We are now ready to establish (|3.3p . Observe that if dist(Q p (x),D) > 0, 
then u XtP = jji f x \U dy and thus (j3.3|) follows from ()A.3P and in this case, 



the constant remains bounded as p approaches r. If u x p = 0, then we have 
dist(Clp(x),D) = 0. In this case, the corkscrew condition (|2,ip implies that 
a(D n n r (z)) > c(r - p) n_1 and (|3~3l) follows from Lemma [Ol 

To establish the inequality (|3.4p we work in a coordinate cylinder and let 
r](x',x n ) be a cut off function which is 1 for \x n — <ft(x')\ < p/2 and for 
\x n — 4>{x')\ > p. We apply the fundamental theorem of calculus and obtain 
that 

(u - u x , p ) a du a _ . x _ drj 

n ^~r^ — \ u - u x, P \ q v+\ u -' u >x,p\ q — 

\u-Ux,p\ oy n dy n 

(A.7) 

We apply the Holder inequality and obtain that the following is an upper bound 
to the right-hand side of t|A.7|) . 



\u— u X)P \ q e n -v da 



I 

Jtt p (x) 



Q p (x) 



du 



dy n 



i/p 



dy 



i/p> 



i(9-iy 



dy 



+C 



Q p (x) 



\u-u x , P r^ r 



If we have the relation (q — l)p' = np/(n — p), or q = p(n — l)/(n — p), we 
may use (|3.3p to obtain 



1/(P'(9-1)) 



It — U 



„n-l 



dy 



< C- 



P) 



n-l 



1/p 



Q r (x) 



\Vu\ p dy 



If we have the relation qn/(n — 1) = np/{n — p), then the estimate (|3.3p will 
give us 

\ 1 /P n -i / » \ ( n - 1 )M 



. In 



u — u 



Qp(x) 



x,p\ 



\qn/(n-l) 



dy 



< C 



(r- p) 



n-l 



Q r (x) 
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Simplifying gives that q and p are related by 1/p = 1/2 + l/(2q) if n = 2 or 
l/p = l/n + (n — l)/(nq) in general, which gives fj3.4f) . 



B The regularity problem in H 1,p (dQ). 

The goal of this appendix is to treat the L p -regularity problem for the Lame 
operator L when the data lies in the Hardy-Sobolev space H l,p (dQ) (see section 
I2.2l for the definition of this space). By the L p -regularity problem, we mean the 
following boundary value problem 

Lu = 0, in n, 

u = f, on dfl (B.l) 

(Vu)* e LP(dfi). 

The boundary values are taken in the sense of non-tangential limits. The argu- 
ment we give is adapted from an argument of Pipher and Verchota [27] used to 
study boundary value problems for the bi-harmonic operator in three dimen- 
sions and was subsequently used to study the traction problem and regularity 
problem for the Lame system by Dahlberg and Kenig |10| for p > 1. Our main 
result is the following theorem which treats the L p -regularity problem in two 
and three dimensions and p < 1. 

Theorem B.2 Let Q C R ra , with n = 2 or 3, be a Lipschitz domain and 
suppose that L is the Lame operator with coefficients satisfying \2. 6\) . 

There exists p\ < 1 so that regularity problem W. 1\) has a solution for 
Pi < P < 1- More precisely we have: 

If f is in H 1,p (d£l), p\ < p < 1, then there exists a solution of the LP- 
regularity problem which satisfies 

du 

\\-Q-\\Hp{dQ) + ||(Vu)*||jy>(an) < C ||/ ||ij-i,p(gn). 

The constant in this estimate depends on M and the global character oftt. 
Furthermore, if p > p\ and u is a solution of the L p -regularity problem, 
IK1\) . with f = 0, then u = 0. 

We recall several facts that will be useful in our proof. A real-variable 
estimate of Dindos and Mitrea |13t Lemma 6.1] tells us that there is a constant 
C so that 

r \ Vp / r \ n/p(n-l) 

/ \u\ p dy) <C( (u*) p(n ~ 1)/n *7 , p>0. (B.3) 

Jn J \Jdtt J 

In work from 1988, Dahlberg, Kenig, and Verchota [llj treated the Dirichlet 
problem for the Lame system (also see work of Mayboroda and Mitrea if n = 2 
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|21j). From this work and a real- variable argument we can show that there 
exists e > so that the L p -Dirichlet problem 



Lu = 0, in 

u = f, on dn (B.4) 

k u* e L p (dn) 

has a unique solution for p in the range (2 — e, 2 + e). This solution satisfies the 
estimate \\u*\\ LP{dn) < C\\f\\ LP{dn) . 

The following technical result is the main step in our argument. The proof 
follows ideas of Pipher and Verchota and especially Dahlberg and Kenig [10} 
Lemma 1.6]. This Lemma will use the sets as defined in section 4. 

Lemma B.5 Suppose f2 is a Lipschitz domain in R", n = 2,3, and L is the 
Lame operator with coefficients satisfying the conditions $2.6}) . If A is a 1-atom 
supported in a boundary ball A r (x), then there exists 7 > so that 

r \ 1/2 

J \Vu\ 2 daj <C2-^ k (2 k r)^ n - 1 ^ 2 , k = 0, 1, 2, . . . . (B.6) 

and for 1 > p > 1/(1 + 7), 

\\(yu)*\\ LPim) <Ca(A r (x))- P -\ 

Proof. We let A be a 1-atom supported in a boundary ball A r (x). From the 
work of Dahlberg, Kenig, and Verchota Theorem 3.7] (see also Mayboroda 
and Mitrea [21] for results in two dimensions), we have a solution of the regu- 
larity problem (jB.ljl with (Vu)* G L 2 (dQ) that satisfies the estimate 



ll(V«)*|| L 2(sn) < C||V t A||^ (Ar(!e)) < Ca(A r (x)y 



-1/2 



This implies the estimate (IB. 6ft with a constant that depends on k. Thus, it 
suffices to prove (|B.6P for k > fco- As in the proof of Lemma f4. 101 we choose ko 
so that dist(A r (x), Sfc ) r. 

To establish (|B.6|) for > ko, fix i? > r and suppose that Agn(y) is a 
boundary ball with Ag^(y) n A r (x) = 0. We use the results Lemma l4~4l and 
(gJU) to obtain 

/ \Vu\ 2 da<^[ \Vu\ 2 dz<-^ [ \u\ 2 dz. (B.7) 

JA R (y) R JO, 2R {y) R Jn 4R (y) 

We let E(y,R) = {z € dfl : T(z) n 4fl (y) / 0} and observe that E(y,R) C 
{z : \y — z\ < CR} where C = C(M, a). We fix an exponent t in (2 — e, 2) such 
that we can solve the L*-Dirichlet problem. The estimate of Dindos and Mitrea 
(IB.3j) . Holder's inequality and the solvability of the L*-Dirichlet problem imply 
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that we have 

(\ (n— l)/n 
W) } J 

< CR n " 3 - ( - n ' 1 

< (7 r ( n - 1 ) 2 /*+ 4 - 2n ^ n - 3 -( n ~ 1 ) 2 /* 



< CR- 3 -^Vf Ltm < C\\A\\l HAr(x)) R^-^ 



(B.8) 



We have used the properties of the atom A to obtain that \A\ < C/r n 2 in the 
last line. If we combine (|R7|) and (|B~8l) we obtain 



_ , n 3-n-(n-l)(2/t-l) 

|V U | 2 do<C(-) o(A R (y))- 1 . (B.9) 

We let 7 = ±((n - 3) + (n - l)(f - 1)). Since t < 2, 7 will be positive if n < 3. 
The estimate (jB.6[) follows from (|B.9[) by an elementary covering argument. 

The second estimate follows using a representation formula for the gradient 
of a solution u as in the proof of Theorem 14.11 Note that it is much easier to 
justify the integration by parts in this case as we have (Vit)* € L 2 (<9£1). I 

Now we are ready to give the proof of the main result of this appendix. 

Proof of Theorem \B.2l Let / = ^ ■ XjAj with each Aj a 1-atom supported in 

Aj = A r .{xj) and J2j ^°'(A?) 1 ~ P < 00 • We let u i be the solution of ([BTT]) with 
/ = Aj as given in the work of Dahlberg, Kenig, and Verchota [11] and consider 
the sum u = Ylj XjVj. Using the estimate (|B.3P and the Poincare inequality, 
we have that as long as p > (n — l)/n, 

-—l 

lbj|| iP n/cn-i)(n) < C\\Vvj\\ L pn/(n-i)^ < C\\{Vvj) ||xj>(an) < Ca(Aj)p . 

The elementary inequality Q2 a j) P — valid if a, > and p < 1, implies 

that for M < M', we have 

A/' / u> N Vp 

II X Villw(n-D(n) < X A i a ( A i) 1_P 

which implies that the series defining u converges at least in L pn /( n ~ 1 )(f2). The 
estimate for ||(V^-)*||iP(gn) hi Lemma fB.51 implies that 

11 cv«)* ib (an) < ^ f E A ^( A i) 1_p ) . p > v(i + 7) 



and that Vit has non-tangential limits a.e. in dfl. We may follow the argument 
in Theorem 16. II to show that du/dp lies in H p (d£l). This completes the proof of 
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existence of solutions to the L p -regularity problem for p > max{(n— 1) /n, 1/(1+ 
7)}- 

We establish the uniqueness of solutions to the L p -regularity problem. Sup- 
pose u is a solution of (IB.ip with u = a.e. on d£l. From Lemma 15,11 or 
the result in Mitrea and Wright [24|, Lemma 11.9] for n = 3, we obtain that 
u* € L q {d£l), l/q = l/(n— 1). Thus, if q > 2-e with 2-e the lower bound 
for uniqueness in the Dirichlet problem from the work of Dahlberg, Kenig, and 
Verchota [IT], we may conclude that u = 0. In two dimensions, this gives 
uniqueness if p > 1 — 1/(3 — e) and in three dimensions we obtain uniqueness 
if p > 1 - e/(4 - e). I 
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